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Abstract 

In spite of the sweeping coherence effects in high energy hadron and nuclei collisions, kj- 
factorization can be recovered for the inclusive gluon production in pA collisions at the leading 
logarithmic order In open charm production at RHIC A:7-factorization holds numerically with 
about 10-20% accuracy. This allows to extrapolate the cold nuclear matter effect observed in 
light and charm meson production in dA collisions to that in AA ones. Unlike the open charm, 
the breakdown of factorization in y/iA production is severe. Indeed, already at the lowest order in 
gluon density the main contribution to the inclusive cross section is proportional to the square of 
gluon density in the nucleus. As a consequence, one cannot infer the cold nuclear matter effect on 
Jjij/ production in AA collisions from that in dA. We present the calculation of J 1 41 multiplicity 
in the framework of the CGC (color glass condensate)/saturation and show that it agrees with the 
experimental data. 



1. Introduction 

One of the cornerstones of the hard perturbative QCD (hpQCD) is factorization theorems, 
which state that soft non-perturbative part of the scattering cross section for any hard process can 
be encoded in universal parton distribution functions (pdf 's). 'Hard' means that the momentum 
transfer Q is much larger than any intrinsic hadron momentum scale. According to Gribov, Levin 
and Ryskin 1 1 1, at high energies there is a universal scale characterizing hadronic wave functions 
- the saturation momentum Qs - that rapidly increases as a power of energy. One therefore ex- 
pects that at higher energies factorization is broken down in a wide region of momenta pj < Qs- 
Understanding the structure of inclusive processes in the region where the hpQCD factorization 
is not applicable is important for quantifying the role of the cold nuclear matter effects in pA and 
AA collisions. In this article I review our present understanding of the subject and its implications 
for inclusive Jj^ production. 

2. Factorization in inclusive gluon and quark production 

Factorization in inclusive gluon production in pA collisions in the low x region was inves- 
tigated in 131, where the cross section was derived that re-sums all leading logarithmic contri- 
butions asln(l/jc) ~ 1 (LLA) for a heavy nucleus in the quasi-classical limit ffjA'^-* ~ 1. One 
does not expect any hpQCD factorization to apply in this case because higher twist interactions 
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of valence quarks and gluons give contributions of order unity. Nevertheless, despite the fact that 
individual diagrams break factorization in covariant and light-cone gauges, the final re-summed 
expression can be cast in the A;7-factorized form. Unlike in hpQCD, the physical quantity that 
is factorized - the unintegrated gluon distribution (f{x, Q^) - is not soft and can be calculated 
perturbatively owing to existence of a hard scale Qs » Aqcd- Another surprising fact is that 
contrary to naive expectations ip(x, Q^) is related not to the momentum space Fourier-image of 
the nucleus gluon field correlation function (A(0) ■ A(x)), but rather to the Fourier-image of 
V^A^(r, b,y), where A^(r,b,3') is the imaginary part of the forward elastic scattering amplitude of 
a color dipole of size r at impact parameter b and rapidity y - ln(l/x) in the heavy nucleus. We 
will see that although the inclusive gluon production in pA collisions is the only known case were 
^/--factorization survives, factorization of the multipoles in the transverse coordinate space is the 
general feature of the low-x cross sections. It must be stressed that this multipole factorization 
does not imply hpQCD factorizations (kj or collinear ones) and neither opposite is generally true. 
A A;/- -factorization formula derived in [2 1 led to successful phenomenology of inclusive hadron 
production in dA collisions at RHIC were the suppression of hadrons at forward rapidities and the 
Cronin enhancement at mid-rapidity were qualitatively predicted ID |4| and then quantitatively 
described in the CGC framework [5 1. 

^T-factorization in pA collisions does not imply any factorization in AA. In fact, no rigorous 
analytic result for an inclusive process exists in the latter case. As in pA, individual diagrams 
in any gauge break the hpQCD factorization. It is however possible that the re-summed result 
exhibit a simpler structure than the individual diagrams. Thus, using a theoretically motivated 
conjecture it has been shown in fT| that inclusive gluon production in AA can be written in a 
form that breaks A:7--factorization only logarithmically. Based on this conjectured approximate 
^7 -factorization, one can derive an important qualitative conclusion about the observed strong 
suppression of inclusive hadron production in AA collisions at RHIC. Since the analogous pro- 
cess in pA collisions in the same kinematic region exhibits Cronin enhancement, we conclude 
that the suppression in AA is not a cold nuclear matter effect. We see that factorization plays a 
key role in establishing existence of a new form of nuclear matter produced in AA collisions. 

Production of heavy quark in pA collisions at low-.*: was calculated in |[8l|9][T0]|llJ. One ex- 
pects that the hpQCD factorization is applicable if the saturation momentum is much smaller than 
the quark mass m 1 12|. At RHIC Qs ~ in for charm and bottom, hence factorization is broken is 
both cases. Indeed, analysis of |,13, | indicates that the semi-classical calculations of ISl fTOl dis- 
agree with A:7-factorization by 10-20% at the f-channel gluon transverse momenta around m. A 
common feature of the inclusive gluon and heavy quark production is that at transverse momenta 
much higher than the saturation momentum, cross sections reduce to the LO hard perturbative 
ones and consequently factorize. In other words, the hpQCD factorization is restored in the kine- 
matic region where the operator product expansion is applicable. It is important therefore that 
the leading term in the twist expansion coincides with that of hpQCD. This is not so in the case 
of the y/t/f production. 



An alternative approach that incorporates the coherence eflfects, but attributes the observed suppression to the pecu- 
liarities of the fragmentation process is presented in |6 1. The two mechanisms can be discriminated by observing whether 
the suppression pattern scales with x^u (CGC) or xo (fragmentation). 



3. Factorization breakdown for inclusive Jlifi production in nuclear collisions 

The mechanism of Jli}/ production in high energy nuclear collisions is different from that 
in hadron-hadron collisions fT4l ITSl |T6l . Consider first the Jjij/ production in pp collisions. 
At high energies formation time of the Jlifj wave function is much larger than the size of the 
interaction region. Indeed, in the nucleus rest frame, the former is 2M^le/{M^' - M^) where 
Ic - 1/{xMn) > Ra is the coherence length and - M^, while the latter is Ra- 

Therefore, we need to take into account only interaction of the cc pair with the nucleus. This 
interaction however depends on the quantum state of the cc pair. In the color singlet model it 
is in the 1 color singlet state while in the color evaporation model it can be in any state with 
invariant mass below the D-meson threshold. Although the color singlet model is physically 
well-motivated, it underestimates the J/t// yield. We therefore assume that cc is produced in the 
1 state with any color Since none of the existing approaches to J /iff production (including the 
non-relativistic QCD model) in pp collisions agrees with all the available data, one may doubt 
the applicability of our model for J/tf> production in pA collisions. Note, however, that there 
are indications that solution to the J/ifr production puzzle lies in understanding the role of the 
higher twist contributions IfTTl . In heavy nuclei, where a^^A^^^ ~ 1, only certain higher twists 
are enhanced - namely those corresponding to interaction of the color dipoles with different 
nucleons. Therefore, higher twist contributions in pA collisions can be taken into account in a 
systematic way. 

Consider now two possible production mechanisms at the lowest order in a, illustrated in 
Fig.[lj (A) g + g ^ J/iff + g, twhich is of the order Oia^A^^^) and (B) g + g + g ^ //lA, which is 
of the order (^(a^A^^^). Since a^A''^ ~ 1, the mechanism in Fig.[IjB is parametrically enhanced. 
Notice, that this leading contribution explicitly breaks A;7-factorization as it is proportional to 

xG(xi)[xG(x2)]\ 




hadron — ha.dro7i collisions hadron — nucleus collisions 



Figure 1: The process of inclusive Jjip production in (A) hadron-liadron and (B) hadron-nucleus collisions. 

At high energies, when Ic » Ra, all possible scatterings of the cc pair inside the nuclear 
medium must be taken into account. The number of scatterings on each diagram is restricted 
however by the requirement that the cc pair is in 1 state. In other words, the number of inelastic 
interactions of the charm dipole in the nucleus must be even. The corresponding scattering 
amplitude reads 

ri„(r,r') = g-'-'ej/s e-'-"G;/8 (cosh[2r ■ r' Q]I%] - l) , (1) 

where r and r' are the cc dipoles in the amplitude and in the complex conjugated one. Of course, 
there is no such restriction on the elastic amplitude involving exchange of gluon pair with the 



vacuum quantum numbers. The elastic amplitude is given by 

rei(r,r') = (l-.-'-^e?/8)(i_,-.-G?/8). (2) 

To obtain the cross section one has two convolute these amplitudes with the virtual gluon and Jij/ 
wave functions, given in e.g. ifTSlfTOll . 

Experimental data is expressed in terms of the nuclear modification factor defined as 
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where iS stands for the overlap area of two colliding nuclei. We calculate the J 1 41 production in 
pp collisions according to Fig.[TJ\ by replacing 7;,, in ([T]i with 



H-l'fQ] 



— e 



(4) 



where no restriction on the number of inelastic interactions has been made. Note also that the 
overall normalization is still not determined since the pp cross section is proportional to the 
probability of soft gluon emission. We approximate this probability by a constant. 

To take into account the low-x quantum evolution we recall that the initial condition for the 
BK II20I l2n evolution equation is given by the Glauber-Mueller formula for the forward dipole- 
nucleus elastic scattering amplitude ll22l 



Nil, b, yo) = 1 - e 



(5) 



Therefore, we can incorporate the evolution effects by writing the scattering amplitudes (|T|i,(j2]),(j4|i 
in terms of the the amplitude and letting the latter depend on rapidity y as dictated by the BK 
equation. Energy evolution of N{i,h,y) is taken into account using the KKT model |5| which 
parameters are fixed to describe the inclusive hadron production in pA collisions. The result of 
the calculation is shown in Fig. [2] We observe that a good agreement with experimental data. 



1.2 F 




Figure 2: Nuclear modification factor for Jjtfj production in iieavy-ion collisions for diff'erent rapidities. Solid (blue) 
line corresponds to RHIC y = 0, dashed (red) lines - RHIC y = 1.7, Dot-dashed (brown) line - LHC mid-rapidity. 
Experimental data from 1231 . 

104 

105 We expect that the breakdown of the hpQCD factorization has important impact on the J I if/ pro- 

106 duction in AA collisions. Because of this the magnitude of the nuclear modification factor Raa 

107 cannot be inferred from the dA calculations. 
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108 4. Conclusions 



109 We discussed the breakdown of the hard perturbative factorization in the gluon saturation re- 

110 gion. This happens due to coherence effects over the entire nuclear length. This effect invalidates 

111 approaches that use pdf's parametrized to include the "nuclear shadowing" corrections. The 

112 very notion of pdf as a universal quantity is not applicable for processes with typical momentum 

113 transfer of the order of the saturation momentum. Fortunately, rigorous analytical perturbative 

114 results can be derived for processes involving scattering off heavy nuclei. In particular, a new 

115 mechanism of Jj^ production in pA collisions has been suggested. It predicts dependence on 

116 nuclear weight A and rapidity y which is in a good agreement with the experimental data. We 

117 plan to extend our calculation to heavy-ion collisions and also calculate the y/i/r polarization. 

118 Our approach can be applied to pp collisions as a model resummation of higher twists. 
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